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ABSTRACT 

We compare QSO emission-line spectra to predictions based on theoretical ionizing continua of 
accretion disks. The observed line intensities do not show the expected trend of higher ionization 
with higher accretion disk temperature as derived from the black hole mass and accretion rate. This 
suggests that, at least for accretion rates close to the Eddington limit, the inner disk does not reach 
temperatures as high as expected from standard disk theory. Modified radial temperature profiles, 
taking account of winds or advection in the inner disk, achieve better agreement with observation. This 
conclusion agrees with an earlier study of QSO continuum colors as a function of disk temperature. 
The emission lines of radio-detected and radio-undetected sources show different trends as a function 
of disk temperature. 

Subject headings: galaxies: active — quasars: general — black hole physics 



1. INTRODUCTION 

The idea of accretion disks around supermassive black 
holes as the power source of active galactic nuclei (AGN) 
is widely accepted. Physical plausibility, the lack of suc- 
cessful alternatives, and the frequent presence of bipo- 
lar jets, support this picture. Qualitatively, the blue- 
ultraviolet "Big Blue Bump" (BBB) in the spectral en- 
ergy distribution of QSOs resembles the expected ther- 
mal emi s sion from the disk atmosphere (Shields 1978; 
iMalkan I 11983(1 . Models achieve some success in fit- 
ting the spe ctral energy distribut ion (SED) of individual 
QSOs (e.g.. iHubenv et all [2001). However, the simple 
disk model has known shortcomings in accounting for 
the soft X-ray co ntinuum, optical polariza tion, variabil- 
ity (see review bv lKoratkar fe Blaeslll999l ). 

In an effort to test t he correspondence b etween disk 
theory and observation. iBonning et ahl (|2007l hereinafter 
B07) examined the optical and ultraviolet colors of QSOs 
in the Sloan Digital Sky Survey as a function of accre- 
tion disk temperature. In simplest terms, do QSOs with 
hotter disks show bluer colors? They found poor agree- 
ment between observed color trends and those predicted 
by standard disk models. For luminosities approaching 
the Eddington limit, the observed colors actually become 
redder with increasing disk temperature, in qualitative 
disagreement with predictions. B07 suggested that the 
discrepancy might result from a modification of the in- 
ner disk within the radius where radiation pressure gives 
a large fractional disk thickness. Modified disk models 
with no emission inside the critical radius gave improved 
agreement with observation. 

Here we extend the results of B07 to include the QSO 
emission lines as diagnostics of the higher frequency por- 
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tion of the disk continuum. The emission lines are gov- 
erned by the intensity of the continuum at ionizing fre- 
quencies. The emitted flux at these frequencies is more 
sensitive to the accretion disk temperature and poten- 
tially offer a stronger test of the trend of the spectral en- 
ergy distribution (SED) of QSOs with disk temperature. 
We again use the large dataset of quasar spectra available 
from the SDSS. We compare the observed quasar spectra 
to model predictions computed with the code AGNSPEC 
(Hubeny et al . 2000) coupled wit h the photoionization 
code Cloudy (|Ferland et al.|[l998l ). With a large range 
of inferred disk temperatures, one might expect strong 
trends in emission-line ratios and equivalent widths as a 
function of disk temperature. We find that these expec- 
tations are not borne out in observed QSO spectra. 

We refer the reader to B07 for further background and 
references to earlier work. We assume a concordance 
cosmology with Ha = 0.7, fi m = 0.3, and H = 70 
km s _1 Mpc _1 . 

2. MODELING AND DATA SAMPLING 

2.1. Accretion Disk Properties 

B07 describe the physical basis of estimating the disk 
effective temperature. In standard thin-disk theory 
(|Shakura fe Sunvaev 1119731 ). accreting matter gradually 
spirals inwards as viscous stresses transport angular mo- 
mentum outward. These stresses lead to local dissipation 
of energy that flows vertically to the surface of the disk 
and is radiated by the photosphere. The resulting effec- 
tive temperature, given by F = a T e g 4 , shows a radial de- 
pendence T e g oc r -3 / 4 at radii substantially larger than 
the inner boundary of the disk. As the inner boundary 
is approached, the effective temperature reaches a maxi- 
mum T max and then drops to zero at the inner boundary 
if the inner boundary condition is one of zero torque. For 
a disk radiating locally as a black body, the disk spec- 
trum is fixed by T m ax, for a given value of black hole 
spin. For realistic opacities, there are significant spectral 
features such as Lyman edges o f hydrogen and heliu m, ei- 
ther in emission or absorption (|Hubenv et al.|[200lD that 
can cause disks with the same value of T max to differ 
in their detailed continuum spectra. However, T max still 
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remains a useful sequencing parameter. 

The key parameters for a QSO accretion disk are the 
black hole mass and spin, the accretion rate, and the ob- 
server viewing angle. An important reference luminosity 
is the Eddington limit, 
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where n e is the electron scattering opacity per unit mass 
and M 8 = M B h/(10 8 M ). The total luminosity pro- 
duced is 



L bo i = eMcf = (10 4t> ' b erg s _1 )e_iM, 
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Here Mo is the accretion rate in Mq yr _1 , and e = 10e_i 
increases from 0.057 for a non-rotating hole to 0.31 for a 
rapidly rotating Kerr hole with angular momentum pa- 
rame ter a* = 0.998 (jNovikov k Thornd [l97l iThornel 
[1971 . For a* = 0.998 
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where 



or alternatively, T max 

L 46 = Lboi/(10 46 ergs -1 ). For a Schwarzschild hole, 
T max is co oler by a fact or of 10 046 for a given Mbh and 
^boi (e.g., Shields 1989). We also define a reference ac- 
cretion rate M Ed = L Ed /c 2 = (10 -0 66 Mq yr _1 )M 8 and 
a dimcnsionless accretion rate m = M/M E( j. Note our 
definition of M E <j in terms of e = 1, so that a value 
m = 3.1 gives L E d for a* = 0.998. 

2.2. Deriving Mbh fl^d T max 

Derivation of Mbh from the width of the H/3 or Mg n 
broad emission lines has become an accepted approxi- 
mation in recent years. The FWHM of the broad lines 
is taken to give the circular velocity of the broad-line 
emitting material (with some geometrical correction fac- 
tor). The radius of the broad-line region (BLR), derived 
from echo mapping studies, incr eases as a function o f the 
continuum luminosity, R oc L r . iBentz et al.l (|2009T ) find 
T = 0.52±0.07, consistent with the value 0.5 suggested 
by ph otoionization physics. Here we use (jShields et al.1 
[20031) 

A/bh = (10 7 - 69 M Q )v 2 3000 L M 5 , (4) 
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where w 30 oo 
L 44 = AL A (5100)/(10 4 

The derived value of T max depends almost entirely on 
the broad line FWHM for a given object. The bolometric 
luminosit y can be est i mated as Lbol = /l X AL,\(5100), 
following iKaspi et al.l (|2000) ) . Using Eq. 0] and this ex- 
pression for Lbol, Eq. |3] becomes 
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We use a bolometric correction fac tor f L = Q 
(|Kaspi et al.|[2000t iRrchards et al. Il2006h . 

2.3. SDSS Spectra and T max Bins 

We use the SDSS DR5 spectroscopic data base for 
"quasars" as our primary observational material. We 
selected AGN using an automated spectrum fitting 
program and imposing quality cuts as described by 
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Note. — Mean properties of the SDSS quasars in the various 
Tmax bins. All quantities are logio, with Mbh hi solar masses and 
M in solar masses per year. Subscript "5100" refers to AGNSPEC 
models with accretion rate adjusted to match observed continuum 
at 5100 A. See text for discussion. 



iSalviander et al.l (|2007l ). These fits give values for quan- 
tities such as the flux and width of the broad H/3 line 
and the continuum flux at various rest wavelengths in- 
cluding A5100. The adopted redshift range is z = 0.1 
to 0.8 so as to keep the [O in] nebular lines out of the 
telluric water vapor region. Approximately one-half of 
the objects fall above the QSO/Scyfert galaxy luminos- 
ity boundary at log j/L„(5100) = 44.44 in units of erg s _1 
. Following B07, we derive Mbh and T max for the indi- 
vidual quasars as described above, and then group them 
into bins with logT max ranging from 5.0 to 5.7 defined 
by windows of 4.95 to 5.05, etc. Table 1 gives the aver- 
age values of log Mbh and log^L t/ (5100) for these bins. 
Also given is the accretion rate and resulting log T max on 
the assumption fi, = 9 and a* = 0.998, based on the bin 
averages for log Mbh and log ^L„(5100). (An alterna- 
tive approach to the accretion rate is discussed below.) 
The values of log Mbh systematically decrease from 8.27 
to 7.20 as logT max increases from 5.0 to 5.7. The mean 
luminosity and accretion rate differs only slightly among 
the bins. As a result, the Eddington ratio increases with 
logTmax from logL/LEd = —1-5 to -0.3 across the bins. 

For each of the bins, we form a composite spectrum 
giving equal weight to all objects based on the frequency- 
integrated flux at earth over the spectrum. These com- 
posite spectra were measured to give values of several 
observational quantities (e.g., the [Ne v] emission-line in- 
tensity) used in the analysis below. Line measurements 
were carried out in IRAI0. Figure Q] gives the compos- 
ite spectra and Figure [2] shows the accretion rate for the 
composite spectra as discussed in the preceding para- 
graph. 

2.4. Accretion Disk Models 

Our goal is to compare basic trends in observed and 
predicted AGN properties as a function of T max . To this 
end, we h ave computed disk m odels using the AGNSPEC 
program (jHubenv e t al. 200l]) in the manner described in 
B07. These models take as input parameters the black 
hole mass Mbh and spin paramter a*, accretion rate M , 
viscosity parameter a, and inclination i of the disk axis 

6 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Fig. 1. — Composite spectra of SDSS quasars binned by T max . The thin vertical line represents the wavelength of the He II A4686 line. 
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Fig. 2. — Accretion rate used in the AGNPSEC models of the 
accretion disk. Lower curve (black dashed) is the accretion rate 
based on an assumed bolometric correction fa = 9. Upper curve 
(red solid) shows the accretion rate required to give the mean ob- 
served A5100 continuum luminosity of the objects in the T max bins. 
For the higher values of T max , the bolometric correction is large 
and the accretion rate and resulting T max is higher than for fa = 9. 
This aggravates the discrepancy between predicted and observed 
emission-line properties. 



to the line of sight. For simplicity, we adopt the com- 
monly assumed value a = 0.1 and f ocus on rapid ly ro- 
tating black holes with a* = 0.998 (|ThorndU97l . We 



take cosi = 0.6 as representative of randomly oriented 
AGN subject to an exclusion of inclinations near edge- 
on because of t he dusty torus postulated in the unified 
model of AGN (jAntonucci fc Milled [l985h . The AGN- 
SPEC code computes the vertical gravity and locally 
emitted flux F(R) in the orbiting frame, and calculates 
for each radius the vertical structure and emergent spec- 
trum with the aid of a pre-computcd grid of models. The 
spectrum observed at infinity is summed over all emit- 
ting radii accounting for all rclativistic effects. In order to 
avoid an excessive number of free parameters, we restrict 
our models to the case a* = 0.998. The spin of the black 
holes in AGN has been much debated, with some au- 
thors arguing for moderate values based on likely merger 
histories or fitting of disk models to the spectral energy 
distributions of individual AGN (e.g.. lCzernv et aLll201lL 
and references therein). 

3. PROBLEMS WITH THE STANDARD DISK MODEL 

Here we discuss tests of the success of the standard 
thin accretion disk model in explaining a variety of ob- 
served properties of AGN. We discuss two different ap- 
proaches to determining the accretion rate of the disk, 
which together with Mbh determines T max . First, we 
follow B07 in assuming a universal bolometric correction 
Jl = 9 as described above. Then we discuss an alterna- 
tive approach (sec Section 3.3) in which M is adjusted in 
the disk models so as to reproduce the observed value of 
^L„(5100) . Figure [2] shows the accretion rate derived for 
the composite spectra in these two different approaches. 
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Fig. 3. — Radial temperature profile for the standard disk model 
and the luminosity fitted model, including the alternative models 
with a modified temperature profile in the inner disk. See text for 
discussion. 



Figure [3] shows the radial temperature profile for these 
models as well as a modified accretion disk model dis- 
cussed in Section 4. 

3.1. Emission-Line Intensities 

As found above, typical values of T max for the accre- 
tion disk (for the parameters of our sample) are of or- 
der 200,000 K, corresponding to a black body peak at 
hv 3fcT rj 1.5/Wh, where hv^ = Ijj is the hydrogen 
ionization potential. Thus, the output of the disk in hy- 
drogen ionizing photons, Q(H°), will be sensitive to the 
value of T max , and the sensitivity will be greater for ions 
with higher ionization potentials. The amount of ionized 
gas and the resulting line luminosities, as well as the rela- 
tive intensities of lines from different levels of ionization, 
should depend strongly on T max . 

3.1.1. H/3 Equivalent Width 

A basic application of this concept is the equivalent 
width of the broad H/3 line. For a photoionized nebula, 
the luminosity o f the H/3 line is given by radia tive recom- 
bination theory (jOsterbrock fc Ferlandll2006f ) as 

L(H/3) = (Q/47r)(eH^/a B )Q(H ). (6) 

Here, f2/47r is the covering fraction for the BLR, e^p 
is the recombination-line emission coefficient, is the 
radiative recombination coefficient for hydrogen to quan- 
tum levels n = 2 and higher, and Q(H ) is the ionizing 
photon luminosity (photons per second). This does not 
allow for collisional excitation or other processes that can 
affect H/3 in the BLR. The H/3 equivalent width is then 
given by EW = L(H^)/L A (4861 A). As T max increases, 
we expect a strong increase in EW(H/3) , because the ion- 
izing frequencies are more sensitive to disk temperature 
than the optical continuum emission. Figure |4] shows the 
predicted EW for the set of AGNSPEC models described 
above, on the assumption of a covering fraction of 0.5 in- 
dependent of T max . The expected increase of EW with 
2max is evident. (The covering fraction of 0.5 was cho- 
sen to give agreement with the observed EW of H/3 at 



Fig. 4. — Predicted and observed equivalent width of broad H/3 
line as a function of T max . The SDSS composites show a systematic 
decrease in EW with increasing T max , in contrast to the strong 
increase predicted by the standard disk model. The modified disk 
model shows a leveling off of the EW for higher T max although 
at an excessively large value; the discrepancy in slope remains at 
lower T max 



log T max = 5.0 and is reasonable in the context of the 
unified model of AGN.) 

Figure [4] also shows the observed trend of the H/3 EW 
as a function of T max for our SDSS composites. The 
observed values show a progressive decrease from from 
120 A to 50 A as logT max increases from 5.0 to 5.7. This 
trend is in qualitative disagreement with the model pre- 
diction. While it is possible that the covering fraction 
changes systematically with T max , an ordcr-of-magnitude 
decrease in f2/4vT would be required from logT max = 5.0 
to 5.7. A search for independent signatures of such a 
trend may be worthwhile. However, it seems likely that 
such an extreme trend in covering fraction would have 
been been noticed in prior investigations. Therefore, the 
discrepancy shown in Figurc[4]is an indication that some- 
thing is wrong with the standard disk model as regards 
the ionizing photon output. This discrepancy resembles 
that found in B07 in which the continuum colors do not 
show the expected trend toward bluer colors at higher 
T max . The systematic decrease in H/3 equivalent width 
with increasing T max gives a new context to the long 
known correlation betwe en the EW and FWHM of the 
broad H/3 emission line (|Boroson fc Greed H992Tl . since 
the derived T max depends mainly on the H/3 broad line 
width (see eq. 5). 

3.1.2. Ionization Ratios 

Typical AGN spectra show emission lines from ions 
with a large range of ionization potentials. In the case 
of ionized nebulae, a useful diagnostic of the ionizing 
continuum is the ratio of the lines of He II and H I. 
In a normal nebular ionization structure, He 11 occupies 
an inner volume of the Stromgren sphere (or a surface 
layer of a slab) whose fraction of the total H + volume 
is proportional to the ionizing photon luminosity ratio, 
Q(He + )/Q(H°). Here Q(Hc + ) is the ionizing photon lu- 
minosity above 4 Rydbergs frequency. The recombina- 
tion line ratio /(He II A4686)//(H/3) is in turn propor- 
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Fig. 5. — Predicted line intensity ratio of He ll/H/3 for accretion 
disk model compared with observed ratio for the narrow lines as 
a function of T max . Note strong underprcdiction of He II for the 
lowest temperatures and the strong trend in the model predictions 
in contrast to the observed intensities. The He II narrow line in- 
tensity is difficult to separate from the broad line profile for higher 
Tinax- See text for discussion. 

tional to the nebular average of AT(He +2 )/iV(H+), giving 
7(A4686)//(H/3) = 1.02 Q(Hc+)/Q(H°). (7) 

This assumes that the nebula is optically thick to the 
Lyman continuum. 

Figure [5] shows that Q(He + )/Q(H + ) does indeed in- 
crease strongly with increasing T max in the AGNSPEC 
models. Figure Q] shows the composite spectra for the 
various bins in logT max . The broad He n line is blended 
with Fe II and difficult to measure. The He II line in the 
narrow line spectrum is well defined for our cooler com- 
posites, in which the broad line widths are large. How- 
ever, it is difficult to separate the narrow from the broad 
components of A4686 for the hotter composites, which 
have relatively narrow broad lines. Nevertheless, two 
points can be be made. A subjective inspection of the 
figure indicates little trend in the narrow He II strength 
as a function of logT max . Also, the observed values of 
narrow He n/H/3 for the lower T max bins are larger than 
would be expected for the predicted disk continuum, ab- 
sent an additional hard component such as a power-law; 
and there is little trend with increasing T max . 

For an observationally more feasible test, we consider 
the narrow line ratio of [Nc v] A3426 to [Nc in] A3869. 
The Ne +4 ion requires even higher ionizing photon ener- 
gies than He +2 and thus the [Ne v]/[Ne in] ratio should 
be even more sensitive to T max than He ll/H/3. However, 
[Ne v]/[Nc in] shows no significant trend with T max (see 
Figure |6]). Thus, this observed line ratio, involving emis- 
sion lines from the NLR rather than the BLR, also fails 
to support the expectation of a harder spectrum in the 
ionizing ultraviolet for higher T max . 

3.2. Soft X-rays 

AGN spectra characteristically include emission in the 
hard and soft X-ray bands. In general, this emission in- 
volves frequencies too high to be the result of ordinary 
thermal continuum emission from a standard accretion 



Fig. 6. — Predicted line intensity ratio of Nc v/Nc in for accretion 
disk model compared with observed ratio as a function of T max . 
Note strong trend in the model prediction that is not reflected in 
the observed intensities. See text for discussion. 

disk. lLaor et all ((19921) and others have noted the in- 
consistency of observed soft X-ray fluxes with the sim- 
ple thin-disk model. A Comptonizing corona, producing 
a power-law spectrum ext ending to high photo n ener- 
gies, is sometimes invoked (jCzernv fe Elvislll987D . How- 
ever, for the higher values of T max that we consider here, 
substantial soft X-ray emission from the inner disk is 
predicted, and there is the possibility of overproducing 
soft X-rays. We have computed T max in the manner de- 
scribed above, for the AGN in the sample described by 
lLaor et all ()1997D . Laor gives X-ray fluxes and values 
of the optical (3000 A) to soft X-ray (0.3 keV) spectral 
index defined by L S /L Q = {v s /v )^ s . Figure [7] shows the 
predicted value of a os as a function of T max for a simple 
set of AGNSPEC models with A/ BH = 10 8 M Q . The in- 
crease in a os with increasing T max is evident. In contrast, 
the observed values of a os in the Laor sample, grouped 
into our same bins in Tmax, show no significant trend 
with T max . In particular, for logTmax > 5.2, the mod- 
els predict a greater value of L S /L Q than observed. This 
is significant, because one might invoke a Comptonized 
component to explain the soft X-rays for low T max ; but 
the hotter disks, in the standard model, over-predict the 
soft X-rays even before any Comptonized component is 
added. While the origin of the X-ray flux of AGN re- 
mains uncertain, it is remarkable that observed values of 
a os show no trend with T max . Thus, we have in a os yet 
another qualitative discrepancy between observed AGN 
properties and the predictions of the simple disk model. 

The dependence of a os on disk effective temperature 
has be en previously di s cussed by other authors. For ex- 
ample, [Derocheilet^n (|2009f ) in their Figure 6 show a os 

versus log[(L bol /L Ed ) 1 / 4 M B H 1 ^ 4 ]- Aside from minor dif- 
ferences in the derivation of Lboi and Mbh, their tem- 
perature axis can be calibrated to our log T max scale by 
adding 7.56. The SDSS objects in their plot span a sim- 
ilar range in T max to ours, and likewise show little sys- 
tematic trend in a os with T max . 

3.3. Models with variable bolometric correction 
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Fig. 7. — Predicted optical to soft X-ray spectral index a oe from 
AGNSPEC mode ls compared with observed values derived from 
ILaor et al.l 11199 71 V Note absence of trend with T max for the ob- 
served values. See text for discussion. 

As noted above and in B07, the use of a universal 
bolomctric correction fa = 9 is not consistent with the 
proposition that disk temperatures may differ consider- 
ably among AGN. An alternative approach is to tune the 
value of M in the AGNSPEC models so as to give the ob- 
served value of ALa(5100) for a given T max bin. We have 
computed a set of models in this fashion, still assuming 
a* = 0.998 and fi = 0.6. We find that for logT max = 5.0, 
the resulting model has a value Lboi giving fa fn 10, so 
that the model differs little from the fa = 9 case. How- 
ever, for the higher T max bins, the value of ,/x implied 
by the AGNSPEC fits increases progressively to a value 
of fa ~ 50 for logT max = 5.7. (Here we use the ter- 
minology T max to refer to the original SDSS quasar bins, 
even though the actual T max of the models is now higher.) 
These models involve very large accretion rates which ex- 
ceed the more conventional M based on fa = 9 by nearly 
an order of magnitude. Such large accretion rates put in- 
creased strain on available accretion supplies. They also 
give Lboi much larger than Lsd, so that the inner radii of 
such disks would likely involve serious modifications due 
to radiation pressure and mass loss, as discussed below. 
Nevertheless, we start our discussion of these models by 
revisiting the several observational tests above using the 
predicted disk spectra from AGNSPEC. 

Because fa increases strongly with T max in these mod- 
els, the actual value of T max exceeds the nominal value 
defining the data bin by an increasing factor as the nom- 
inal logTmax increases. These models are useful, never- 
theless, because we may think of the nominal logT max 
spectrum bins as a sequencing procedure that gives us 
a set of SDSS AGN involving a large range in T max , 
even if there are ambiguities as to how well the nomi- 
nal value agrees with the actual disk tcmperaurc. Be- 
cause the models adjusted to T5100 have a greater range 
of M and logT max than the nominal models, they are 
expected to show even stronger trends in observational 
quantities sensitive to T max , such as continuum color and 
high ionization emission-line strength. 

Figurc|4]shows the predicted equivalent width of H/3 for 



the luminosity-fitted models. The strong increase with 
the bin value of T max is now even more pronounced and 
in greater contrast to the observed decrease. Figures [5] 
and [6] show the predicted He 11 and [Ne v] intensity. The 
strong increase in [Ne v] with T max is now in even greater 
contrast to the observed constancy of [Ne v]/[Ne ill]. 
Figure [7] shows the extreme trend in a os versus T max for 
the luminosity-fitted models. The hotter cases greatly 
overpredict the observed soft X-ray fluxes. 

Clearly, the luminosity-fitting approach only exacer- 
bates the discrepancy between the theoretical predic- 
tions and observations. However, this approach may 
have merit physically, since it avoids use of a univer- 
sal bolometric correction. Moreover, for the higher val- 
ues of T max , the fa = 9 models severely undcrpredict 
the luminosity at A5100. If the accretion rate is not 
increased, some restructuring of the AGN continuum 
source is needed to channel energy more efficiently into 
the continuum at optical wavelengths. In any case, the 
standard thin-disk model requires some kind of modifi- 
cation to explain the continuum and emission-line obser- 
vations. 

4. A MODIFIED DISK MODEL 

B07 found that observed continuum colors of SDSS 
QSOs did not become progressively bluer with increasing 
Tmax in the manner predicted by the AGNSPEC models. 
We suggested that this resulted from higher Eddington 
ratios at higher T max , leading to departures from the 
standard thin disk model. This might involve thickening 
by radiation pressure, inefficient "slim disk" accretion, 
and mass loss in a radiation-driven wind (see also De- 
roches et al. 2009). As an estimate of the potential effect 
on the continuum colors, we computed AGNSPEC mod- 
els in which the inner disk was truncated inside the radius 
where radiation pressure gave a vertical scale height H 
for the disk, relative to the radius, of H/R = 0.5. This 
led to improved agreement with the observed colors, in 
particular giving a non-monotonic trend of color with 

J max • 

Here we explore a more physical model in which the 
disk effective temperature is modified inside the criti- 
cal radius. The basic idea, following IPoutanen et all 
(120071 ) and lBegelman. King, fc Pringlej ([20061 ) . is to take 
account of the luminosity produced at a given radius rel- 
ative to I/Ed- If the luminosity emitted by a wide an- 
nulus, dlnL/dlni?, is well below Lsd, then the disk will 
not be severely thickened by radiation pressure, and stan- 
dard disk physics applies. However, except near the inner 
boundary, dlnL/dlni? increases with decreasing R. For 

sufficiently large values of M/Mbh, a radius is reached 
where dlnL/dlni? reaches LEd- Near this radius, serious 
modifications of the disk s tructure can be anticipated. 
IShakura fc Sunvaev I (|1973f ) called this the "spherization 
radius" R sp on the basis t hat radiation pre s sure w ould 
severely thicken the disk. IPoutanen et al.l (|2007| ) find 
that, inside i? sp , a combination of wind mass loss and 
advection occurs to a degree that keeps the local lumi- 
nosity dlnL/dlni? close to Lsd- This results in a break 
in the run of effective temperature with radius, so that 
inside i? sp one has T e g cx i? -1 / 2 rather than the stan- 
dard dependence T e g cx i? -3 / 4 . This reduces the value 
of Tmax actually reached in the inner disk as well as the 
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luminosity radiated from the hottest parts of the disk. 
The value of R sp in units of the gravitational radius is 

approximately r sp = R sp / R g ~ rh. As M increases, the 
effective temperature at R sp actually decreases, following 



sp 



(10 b 



(8) 



Here we have used the Newtonian expression for the lo- 
cal flux, omitting the normal factor of 3 increase in flux 
due to viscous transport of ene rgy from smalle r radii . 
This modification was found by IPoutanen et al.l (j2007T ) 
because of energy loss in the inner disk due to advection 
and mass loss. For large values of m, a luminosity ~ L Ed 
is emitted in the vicinity of i? S p, and a similar luminos- 
ity is emitted per unit In R at smaller radii. Thus T sp 
substantially controls the character of the emitted spec- 
trum. Hence, the fact that T sp decreases with increasing 
rh once spherization comes into play is suggestive of the 
reversal in the trend of colors with T max found in B07. 
(Note that, in this context, Tmax is defined in terms of 
the accretion rate entering the disk at large radii, and dif- 
fers from the actual maximum temperature in the inner 
disk.) 

The potential for an inverse relationship between M 
and disk color temperature is strengthened if one con- 
siders the optical depth of the disk wi nd. In the case 
of hig hly super-Eddington accretion, IPoutanen et al.l 
(|2007 ) find that a large part of the inflowing material 
is expelled as a wind in the vicinity of R sp . In the simple 
approximation of a spherical wind with constant radial 
velocity equal to the circular orbital velocity at i? S p, the 
optical depth of the wind, measured radially from R sp , 
is r. 



,V2 . 



The optical depth in the wind, mea- 
sured outward from some radius R that exceeds R sn , is 



(> R) 



For the radius of the photo- 



sphere, where r(> R) = 1, this gives r p h ~ rh 3 / 2 . 
effective temperature at the photosphere is 



P h 



(10 5 - 75 K)m 



- 3 / 4 m 8 - 1/4 . 



The 



(9) 



This raises the prospect that, for super-Eddington ac- 
cretion rates, the maximum visible disk temperature will 
be T p h rather than T max . An estimate of the progres- 
sion of the disk temperature as it actually appears to the 
observer is then 



T app — min(T max , T p h) 



(10) 



This assumes that there is enough absorption opacity in 
the wind to reprocess the emerging energy; otherwise, 
the escaping spectrum may still be characterized by T sp . 

For our Series 1 models above, with an accretion rate 
based on fz = 9, the value of logL/L^d reaches -0.3 
for logT max = 5.7. Thus, unless opacity sources other 
than electron scattering reduce the effective Eddington 
limit, the effects of radiation pressure may be modest. 
However, for our Series 2 models with M adjusted to 
fit the continuum at A5100, logL/LEd reaches +1.1 
for logT max = 5.7. For this series, at logT max = 
(5.0,5.4,5.7) we have logm = (-0.79, +0.29, +1.65); 
logmin(T max ,T sp ) (5.05,5.44,5.13), and 

logmin(T max ,T ph ) = (5.05,5.37,4.75). Thus, a re- 
versal of the trend of the apparent temperature of the 
disk as a function of the nominal T max would be expected 
at logT max 5.4. The fact that this value agrees with 



the turn-around in the color trend found in B07 is 
fortuitous, given the approximations in our discussion. 
However, this simple calculation does illustrate the 
possibility of reverse trend of color temperature with 
M/Mbh 2 for super-Eddington rates. 

In order to explore the impact of super-Eddington 
accretion on the disk spectrum and resulting emission- 
line intensities, we have computed modified models with 
AGNSPEC in which the standard T g inside a break ra- 
dius Rb is multiplied by a factor (i?/i?f,) 1/4 (see above). 
For the value of the break radius, we took Rb/r g = 
Xb (L/Lea). We present here results for Xb = 7, which 
was motivated by a preliminary analysis but which serves 
here simply as an illustrative example. We used the same 
values of Mbh and M as used in our Series 1 and Series 
2 models above. The run of T e ff with radius for the mod- 
ified models is shown in Figure [3] As expected, there 
is little effect for Series 1 (fz = 9); but substantial al- 
terations are seen for Series 2. In Figure IU the Series 2 
prediction for the H/3 equivalent width levels off above 
logT max = 5.4. This reduces but does not eliminate the 
discrepancy with the observed trend. Figure [5] shows a 
similar leveling of the predicted He II intensity. Figure 
[6] shows a substantial reduction in the predicted Nc v 
intensity for logT max around 5.4. These models are not 
intended to achieve a detailed fit to the observed spectra, 
but they do illustrate the substantial effect on the ioniz- 
ing continuum that results from a modified disk structure 
resulting from super-Eddington accretion. 

5. RADIO VS. NON-RADIO SOURCES 

Our focus has been on the entire SDSS quasar data 
set, which is strongly dominated by radio quiet objects. 
However, radio loudness raises interesting questions re- 
garding the black hole spin, which affects the inner disk 
boundary and thus the maximum temperature reached 
in the disk. Numerous authors have discussed the possi- 
bility that r adio loud AGN h ave rapidly spinning black 
holes (e.g.. iBlandfordl I1990T ) . and on the other hand 
IGarofalo et al.l (12010ft have proposed that radio loud ob- 
jects have retrograde black hole spin with respect to the 
disk angular momentum. The innermost stable circular 
orbit (ISCO) decreases from QR g for a = to 1.22R g for 
a* = 0.998. Correspondingly, the radius of maximum ef- 
fective temperature moves inward from 9.5R g to 1.55i? 9 ; 



and the value of T max increases from 10 
.5.43 7|>l/4 iu-i/4 



4 - 79 M 1/4 M 8 - 1/4 



to 10 5 - 43 M V M s 1 (|Shieldslll989l) . For retrograde spin, 
the ISCO moves still farther out. These differences in the 
inner disk temperature could have a significant effect on 
the ionizing spectrum of the disk, while having less ef- 
fect on the optical luminosity. Thus, we might expect 
differences between radio loud and quiet AGN in the ob- 
servational quantities sensitive to T max discussed above. 

We have formed composite spectra separately for the 
radio loud and radio quiet AGN, adhering to our previ- 
ously defined bins in T max . Here radio loud was taken 
simply to mean d etected in the FIRST radio survey 
(|Becker et al.lll995[). This is not equiva lent to the stan- 
dard definition (jKellermann et al.lll989| ) and is distance- 
dependent, but it is sufficient to look for any qualita- 
tive difference between radio loud and quiet AGN. Fig- 
ure [8] shows the observed trend of the [Ne v] intensity 
with T max for the radio loud and quiet composites. The 
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Bonning et al. 
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data in terms of T max helps to clarify what trends are ex- 
pected. We find that the predicted trends in the intensity 
of H/3, He II and ]Nc v] arc not seen in the observations. 
Over a factor of five in the value of T majX as inferred from 
Mbh and fL„(5100) , the [Ne v] intensity shows no sys- 
tematic variation, and the H/3 equivalent width shows a 
reverse trend. While our models are highly simplified, in 
particular for cases that approach or exceed the Edding- 
ton limit, the basic immunity of these key line ratios as 
well as a os to the expected value of T max is a striking 
empirical fact. These results suggest that a serious re- 
examination of the structure and emission mechanisms 
of the AGN central engine is in order. This might rea- 
sonably start with more detailed studies of the nature 
of disks in systems with high Eddington ratios, taking 
account of advection, winds, and all contributing sources 
of opacity. 



Fig. 8. — Comparison of the [Ne v]/[Nc III] ratio for radio de- 
tected and non-detected SDSS quasars as a function of T max . Note 
significant difference for the higher values of T ma x- See text for dis- 
cussion. 



[Ne v] /[Ne m] ratio is similar for the two groups for 
logT max = 5.0, but there is a decline for the radio de- 
tected objects with increasing T max . Certainly, there is 
no evidence that [Nc v] is stronger for the radio-detected 
objects, as might be expected if they have higher spin of 
a prograde sense. On the other hand, the weaker [Ne v] 
of the radio objects might be expected if they involve 
retrograde disks. 

6. CONCLUSION 

The results presented here give focus to concerns about 
the disk model for the AGN continuum that have been 
voiced in the literature by numerous authors. In general 
terms, AGN over a wide rang e of Mrh and lumin osity 
have rather similar SEDs (see lLaor fc Davisl 1201 lL and 
references therein). Organization of the observational 



We thank Ivan Hubcny for the use of the AGN- 
SPEC program, R. Antonucci, O. Blacs, A. Laor, and 
B. Wills for helpful discussions. G.S. gratefully ac- 
knowledges the support of the Jane and Roland Blum- 
berg Centennial Professorship in Astronomy. Fund- 
ing for the Sloan Digital Sky Survey (SDSS) has been 
provided by the Alfred P. Sloan Foundation, the Par- 
ticipating Institutions, the National Aeronautics and 
Space Administration, the National Science Foundation, 
the U.S. Department of Energy, the Japanese Monbuk- 
agakusho, and the Max Planck Society. The SDSS 
Web site is |http: / /www.sdss.org/| The SDSS is man- 
aged by the Astrophysical Research Consortium (ARC) 
for the Participating Institutions. The Participating 
Institutions are The University of Chicago, Fermilab, 
the Institute for Advanced Study, the Japan Participa- 
tion Group, The Johns Hopkins University, the Korean 
Scientist Group, Los Alamos National Laboratory, the 
Max-Planck- Institute for Astronomy (MPIA), the Max- 
Planck- Institute for Astrophysics (MPA), New Mexico 
State University, University of Pittsburgh, University 
of Portsmouth, Princeton University, the United States 
Naval Observatory, and the University of Washington. 



REFERENCES 



Antonucci, R. R. J.; Miller, J. S. 1985, ApJ, 297, 621 
Becker, R.H., White, R.L., & Helfand, D.J. 1995, ApJ, 450, 559 
Begelman, M. C, King, A. R, & Pringle, J. E. 2006, MNRAS, 
370, 399 

Bentz, M. C, Peterson, B. M., Netzer, H., Pogge, R. W., & 

Vestergaard, M., 2009, ApJ, 697, 160 
Blandford, R. D. 1990, in Courvoisser, T. J.-L., Mayor M., eds., 

Active Galactic Nuclei. Springer, Berlin, p. 161 
Bonning, E. W., Cheng, L., Shields, G. A., Salviander, S., & 

Gebhardt. K. 2007, ApJ, 659, 211 
Boroson, T. A., & Green, R. F. 1992, ApJS, 80, 109 
Czerny, B., & Elivis, M. 1987, ApJ, 321, 305 

Czerny, B., Hryniewicz, K., Nikolajuk, M., & Sadowski, A., 2011, 

MNRAS, 415, 2941 
Deroches, l.-B., Greene, J. E., & Ho, L. C. 2009, ApJ, 698, 1515 
Ferland, G. J., Korista, K.T., Verner, D.A., Ferguson, J.W., 

Kingdon, J.B. & Verner, E.M. 1998, PASP, 110, 761 
Garofalo, D., Evans, D. A., & Sambruna, R. M. 2010, MNRAS, 

406, 975 

Hubeny, I., Blaes, O., Krolik, J. H., Agol, E. 2001, ApJ, 559, 68C 
Kaspi, S., Smith, P. S., Netzer, H., Maoz, D., Jannuzi, B. T., & 

Giveon, U. 2000, ApJ, 533, 631 
Kaspi, S., Maoz, D., Netzer, H., Peterson, B. M., Vestergaard, 

M., & Jannuzi, B. T. 2005, ApJ, 629, 61 



Kcllcrmann, K. I., Sramek, R., Schmidt, M. Shaffer, D. B., & 

Green, R. F. 1989, AJ, 98,1195 
Koratkar, A., & Blaes, O. 1999, PASP, 111, 1 
Laor, A., & Davis, S. 2011, preprint arXiv:1110.0653 
Laor, A., Fiore, F. Elvis, M., Wilkes, B. J., McDowell, J. C. 1997, 

ApJ, 477, 93 
Malkan, M. 1983, ApJ, 268, 582 

Novikov, I. E., & Thorne, K. S. 1973, in Black Holes, Lcs 

Houches, ed. C. DeWitt & B. DeWitt (New York: Gordon & 
Breach) . 

Osterbrock, D. E. & Ferland, G. J. 2006, Astrophysics of 

Gaseous Nebulae and Active Galactic Nuclei (2nd ed.; 

Sausalito, CA: University Science Books) 
Poutanen, J., Lipunova, G., Fabrika, S., Butkevich, A. G., & 

Abolmasov, P. 2007, MNRAS, 377, 1187 
Richards, G. T., Lacy, M., Storrie-Lombardi, L. J., Hall, P. B., 

Gallagher, S. C, Hines, D. C, Fan, X., Papovich, C. et al. 

2006, ApJS, 166, 470 

Salviander, S., Shields, G. A., Gebhardt, K., & Bonning, E. W. 

2007, ApJ, 662, 131 

Shakura, N. I., & Sunyaev, R. A. 1973 A&A, 24, 337 
Shields, G. A. 1978, Nature, 272, 706 

Shields, G. A. 1989, New York Academy Sciences Annals, 571, 110 
Shields, G. A. 2003, et al., ApJ, 583, 124 



Accretion Disk Temperatures of QSOs 

Thome, K. S. 1974, ApJ, 191, 507 



